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1. Introduction 



The purpose of this paper is to give a mathematical explanation of a formula for 
the scattering matrix for a manifold with infinite cylindrical ends or a waveguide. 
This formula, which is well known in the physics literature, is sometimes referred to 
as the Mahaux-Weidenmuller formula [9]. We show that a version of this formula 
given in (jl.7p below gives the standard scattering matrix used in the mathematics 
literature. We also show that the finite rank approximation of the interaction matrix 
gives an approximation of the scattering matrix with errors inversely proportional to 
a fixed dimension-dependent power of the rank. 

Theorem 1. Let X = Xq U (0, oo) x OXq be a manifold with cylindrical ends - 
see §m/br a precise definition and Figured for an illustration. Let {ty n }™ =0 be an 
orthonormal set of real eig en functions of the Neumann Laplacian, —H- m , on X with 
eigenvalues —t 2 . Let {<£>a} be the same set for the Laplacian on dX , with —a 2 
denoting the corresponding eigenvalues. Let us define the interaction matrix by 

W N , A {k) : L 2 (dX ) — >L\X ), 
( L1 ) W N , A f= E (&-°l) i (*n\ox ,<Px){f,<Px), 

0<t„<VN 0<ctx<VA 

and the effective Hamiltonian by 

H N>A (k) = H in - iW N , A (k)W N , A (ky . 

Then for fcGl, the entries of the scattering matrix (see £TJ|) are given by S\ t y(k) = 

(1.2) - ((/ - 2iW N , A {k)\k 2 - H^k^W^ik))^,, <^ A > + 0(JV-s + e~ k / c ) , 

if a\,a\/ < \k\, and A > k 2 . The error bound 0(N~^) is optimal - see £0 and the 
constant can be chosen uniformly for k lying in compact sets. 

Theorem [2] provides a related result, for other values of k. Also, we remark that the 
matrix defined by the leading term in (jl.2p . o~\,ay < \k\, is in fact unitary - see 
Lemma 15.11 

l 
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The physics literature contains several versions of the Mahaux-Weidenmiiller for- 
mula. One commonly found formula - see for instance p],[ll] and references given 
there - is given as follows 

(1.3) S f (k) = - (i - 2iW{k)*{k 2 - H^Wikfj . 
Here 

(1.4) H cS d ^ H m - iW(k)W(k)* 

where — Hi a is the Neumann Laplacian in the "interaction region" Xq, a compact 
piece of the waveguide or manifold with infinite cylindrical end, and W(k) is the 
frequency dependent interaction matrix. When applied in numerical simulations only 
finite number of modes of H in are taken which results in a finite rank approximation 
of W(k), as described in (11.11) . The formula, in its finite rank version, is the basis of 
random matrix models in scattering theory - see [6j Section III.D]. For some recent 
experimental results related to the formula see for instance |14j . 

The formula (jl.3p is not strictly speaking correct. The advantage of (11.31) is that 
Sf(k) is unitary for real k by a linear algebra argument. It is also close to the correct 
scattering matrix given below. 

As shown in Proposition I3.5[ the scattering matrix [2] which is standard in the 
mathematical literature is recovered from an expression close to (11. 3ft : 

(1.5) S f (k) = - (I - 2iW{k)\k 2 - H ca y l W(k)) 
with 

(1.6) H cS = H in -iW(k)W(k)\ 
and, with the notation of (11.11) . 

W(k) = f Woo :00 (k) . 

In fact, —H e ff = —H c ff(k) is the Laplacian on Xq, with a boundary condition that 
depends on fc; see Lemma 13.21 Lemma 13.31 demonstrates the relationship between 
(k 2 — iJ c fj) _1 and the resolvent of the Laplacian on X. 

This correct version (11.51) appears in p], though again only a finite number of 
modes are included. We note that our sign convention, while agreeing with [1], is 
not consistent with many other authors. It appears that this sign is correct, and 
that the difference can be traced to a different normalization of the scattering matrix. 
The difference between fl 1 . 3 [) and (11.51) does not appear in many of the physics papers, 
where generally only an approximation W a (k) of W(k) is used, and the approximation 
is such that W a {k)* = W a (ky. The operator Sf(k), unlike Sf(k), is typically not 
unitary for real k. 
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However, (11.51) gives what one might call the extended, or full, scattering matrix. To 
get the usual finite dimensional unitary scattering matrix (whose dimension changes 
at roots of the eigenvalues of the cross section of the end), we put, for k real, 



most k 2 . Here Ay is the Laplacian on the cross section of the end. Proposition 13.51 
shows that this is the unitary scattering matrix which appears in the mathematical 
literature. Lemma 15.11 gives an algebraic proof that the matrix given by (11.71) is 
unitary for k G JR.. Note that if k £ R, the operator defined by ( 11.31) is unitary, but 
the finite rank-operator (corresponding to a finite-dimensional matrix) 



with H e ff given by (ll.4p . is not unitary in general, if W(k) takes into account contri- 
butions of evanescent modes. Evanescent modes correspond to eigenvalues of —Ay 
larger than k 2 . 

Let us add that the articles [1] and [11] already have a fairly mathematically care- 
ful description of the Mahaux-Weidenmuller formula. In [12] a detailed analysis of 
several one dimensional models is also provided. Another related approach to scat- 
tering/transport is due to Fisher-Lee [5], see also [3]. 
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aging us to write this paper, Henning Schomerus for letting us know about the Fisher- 
Lee formalism, Ulrich Kuhl for helpful conversations, and an anonymous referee whose 
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its hospitality in spring 2009. 

Remark. We use the notation (u, v) to denote the Hermitian inner product, and 
(u, v) to denote the form which is linear in both arguments. 



In this section we recall the general assumptions for manifolds with cylindrical ends 
and the definition of the scattering matrix. 

Our model is a manifold X with infinite cylindrical ends and smooth metric g- see 
Figure [TJ In physics language that means a waveguide with periodic boundary con- 
ditions. The same arguments apply to waveguides with Dirichlet or Robin boundary 
condition but we choose to avoid mild technical complications associated with that 





2. Scattering matrix 
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X 



Figure 1. An example of a manifold with an infinite cylindrical end. 



setting. For purely notational reasons we also assume that there is only one end. 
Then 

X = X U (0, oo) x Y , Y = dX , 

where X is a compact manifold with a smooth boundary Y. We require that 
9 \[o,oo)xy= (dx) 2 + gy, where gy is a metric on Y. Moreover, we choose our de- 
composition so that there is a neighborhood U C X of dX on which g is a also a 
product: 

g\u= (dx) 2 + g Y . 

Recall that {<£>a} are an orthonormal set of eigenfunctions of Ay. We use the 
convention that the energy is k 2 , and k\ = a/A; 2 — cr|, with the imaginary part chosen 
to be non-negative when ImA; > 0. We call the region with ImA; > the physical 
region. Given A G N, if k is in the physical region, and with Im/c > 0, there is a 
unique $a(p, k) so that 

(-A x -fc 2 )$ A (p,fc)=0onX 



(2.1) 
and 

(2.2) 



r ik xX fM + J2S yx {k)e 



Ik 



for some Sy\. To see this we use the resolvent (-A x — k 2 ) 1 which is a bounded 
operator L 2 (X) -> H 2 (X), for ImA; > 0: 

$ A (p, k) = (1 - i^A(l/)e- tt *" + (-Ax - fc 2 )' 1 ([A x , ^] (^ A ( 2 /)e-^)) (p) , 

Since on Xi we have —Ax = —d 2 — Ay, separation of variables shows that $a can 
be written as in (12.21) . 
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The resolvent, (—Ax — k 2 ) 1 , continues meromorphically to 
(2.3) K(A 9Xo )D{k : lmfc>0}, 

a Riemann surface branched at <ta's - see [TUl Sect. 6. 7]. We remark that this Riemann 
surface is such that each k\ defined above extends to be a holomorphic single-valued 
function. Thus &\{p, k) has a meromorphic continuation to A a ^ gx ) which is regular 
for Irak = except when fcy's are 0, or when k 2 G u(— Aj). 
The full, or extended, scattering matrix is the infinite matrix 

Sf(k) = (Sx>\(k))x,\>eN- 

For k G K, the matrix more commonly called the scattering matrix is the finite- 
dimensional matrix given by 



S(k) = (Sx'x(k)) a 2 



,<k 2 - 



We remark that if Imk > 0, while each entry Sxy{k) is well-defined away from its 
poles, there is not a canonical choice for "the" scattering matrix. However, in general 
it is (\^/y/ky)Sx'x(k), not SVa, which has a meromorphic continuation to A (T (A aA - o ) 
for each A, A'. We shall use this continuation in the proof of the theorem. 

3. The formula 

Let Ay be the Laplacian on Y, and let {cr\} be the eigenvalues of —Ay, repeated 
according to multiplicity, and let {v?a} be an associated set of real, orthonormal 
eigenfunctions of the Laplacian on Y . Let — H in be the Laplacian with Neumann 
boundary conditions on Xq, and let {\l/ n } be a set of real, orthonormal eigenfunctions 
of H in . 

First, we define the operator W(k) by explicitly giving its Schwartz kernel. Our 
starting point is the representation of W(k) from [TJ or [TTJ. We write p to represent a 
point in X , and y or y' to represent a point in Y; on U C X we may write p = (x, y), 
with {x = 0} = dXo- Then, with 

*„,A(o) d =y <f X (y)Mo,y), 

we follow the physics literature and define the coupling operator by giving its integral 
kernel (with integration with respect to Riemannian densities) as 

w( P ,y') = Vh^niP^uMMy') 

(3.1) n ' X 

= J2^n(p)Pk^n(0,y'). 

n 

Here = (k 2 + Ay) 1//4 is defined by Pk^Px = V^a^a- While either choice of the 
square root is possible, it is crucial that this is consistent with that used to define 
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the scattering matrix; see (I2.2p . The series converges in the sense of distributions. 
Hence, W(p,y') is understood as a distribution on X x Y - see Lemma [3.11 below. 

This definition (13.11) appears normally in the physics literature. When we take all 
the eigenstates then W (and especially W l ) takes a very simple form given in the 
following lemma. When, as is done in the physics literature, we take only finitely 
many states, the formula for W* is given in the remark after the lemma. 

Let V'(X ) be the space of distributions on X supported in the (closed) set X . 
We also use the following convention: for s > the space H S (X ) denotes restrictions 
of elements of H S (X) to X Q , while for s < 0, H S (X ) denotes elements of H S (X) 
supported in the (closed) set X Q . See Appendix B.2] for a careful discussion: in 
the notation used there 

H {Xo) ~ { H (S) (X ) s<0. 
With this notation in place we can formulate 

Lemma 3.1. The operator (13. ip is equal to 

(3.2) W(k)g = 6 9Xo P k g, g G C°°(dXo) , 
where 5sx G ^'(^o) is the distribution defined by 

$dx (p) = / <f\dx rfvoly . 

JdX 

We have 

(3.3) W(k) : H S (Y) -> J ff min (" 1 / 2 "' S " 1 )(X ) , sel, 

and W(k)g\x S = 0. The transpose, W(ky : H x / 2+s {X Q ) -> H~ 1 / 2+s (Y) ! s > 0, is 
given by 

(3-4) W(kYf(y) = P k (f\ dXo ). 

Proof. To prove (13.21) we need to compute, in the notation of distributions, W(f ®g), 
where / G Cq°(Xq). The definition (13. ip gives 

W(f ® g) = V ( I * n f) ( [ Vnhxo P k g 

„ \JXo J \JdX 



7./ * n r 9 x 1 p*» 

/tax 

'8X 

which proves (13. 2p and, by duality, (13.41) . The mapping property of W{k) 1 follows 
from the fact that f ^ f \ Xo takes H S+1 / 2 (X ) to H s (dX ) for s > 0, and P k : 
H s {dX ) -> H s - l / 2 (dXo). The mapping property Q follows by duality. □ 
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Remark. In Lemma 13.11 all the structure of the the basis of eigenvectors of H in 
and Ay disappears. The question which we address in Section H] is how close the 
approximation based on using only finitely many basis elements gets to the actual 
scattering matrix. Then for a = (N, A) G [0, oo] 2 we define 

W a (kY = P k l [0iN] (A Y )R i[0,A](#in) , RU = U\ 8X0 . 

We note that 

W(oo,oo)(&) = W(k) , 
and that for N < oo and A < oo, 

W a (k) : V{Y) — C°°(X ) , W a (kY : V'(X ) — > C°°{Y) , 

where C°°(Xq) denotes extendable smooth functions on the compact manifold X . 

We make the definition (11.61) of H c r rigorous via the quadratic form 

q{u,v) = q(k)(u,v) = / VuVv dvo\ Xo - i \ W t (k)uW*(k)v dvo\ Y 

JXo JdX 



VuVv dvo\ Xo — i / PkRuP^Rv dvoly 
x JdX 

with form domain H 1 (X ). If 

q(u, v) = (w, v) 

for some w G L 2 (X ) and all v G H 1 (Xq), then u is in the domain of H c g and 
H c ffU = w. Moreover, 



(w, v) = q(u, v) = — / A Xo uv + / (d n u — iP k Ru)v , 
Jx JdX 

where d n u denotes the outward unit normal derivative at the boundary. Since this 
must hold for all v G H l (X ), —A Xo u = w and 

= d n u - iPlRu. 

We note that u G H 2 (X ) where the space is defined by restricting elements of H 2 (X) 
to Xq - see [71 Appendix B]. We summarize this in the following 

Lemma 3.2. Suppose u G Domain(iJ c fj). Then u G H 2 (X ), and 

H e ffU = —A Xo u , d n u — iP^Ru = . 

Next we investigate the relation between (k 2 — if c ff) _1 and the resolvent of the 
Laplacian on X. Denote 

R x (k) = {k 2 + Ax)" 1 , for ImJfc > 0. 

Then, for K C X any compact set \xRx{k)\K has a meromorphic extension to 
A (T (a S x ) ; see M- In LemmalOwe shall show that (k 2 - H^y 1 : L 2 (X ) H 2 (X ) 
exists for k 2 0, Irak > 0, and is meromorphic on A CT (A ax ), the Riemann surface 
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(12. 3p . One could provide an alternate proof using the first part of the proof of Lemma 
13.31 and the results of [10] on the meromorphic continuation of Rx{k). 

We remark that when we use k G A ct (a SXo ), by abuse of notation we mean by k 2 
the complex number which is the continuation of k 2 from the physical half plane 
Im k > 0. 

Lemma 3.3. We have the following relation between Rx{k) as defined above and 
(k 2 - H eS )^ = (k 2 - H eS (k))- 1 : 

(k 2 -Hcs)" 1 = l Xo R x (k)t Xo 

for k G A CT (A aXo ). In particular, the poles of (k 2 — if e fj) _1 are the same as the poles of 
Rx(k). 

Proof. Suppose g G L 2 (X ) C L 2 (X) and g is in a neighborhood of dX . Then 

(3.5) (k 2 + A Xo )l Xo R x (k)g = (k 2 + A Xo )l Xo R x (k)t Xo g = g on X° 
for k G K(A eXo )- 

Note that since supp g C Xq, (k 2 + A x )g = on X±. Then for ImA; > (that is, 
for k in the physical space), the requirement that R x {k)g G L 2 (X) means that 

(3.6) R x (k)g\ Xl =J2 a * eikXXl P\ 

for some constants a\ = a\(k). But then, using the support conditions of g there is 
a neighborhood U C X of dX so that 

Rx(k)g\u=J2 a ^ kxX ^- 

Thus 

{d n -iP 2 k )(R x (k)g\ Xo )\ dx = Q 

so that R x {k)g\ Xo is in the domain of H cS = H cS (k). Together with (I3.5p . this means 
that 

(k 2 -H cS y 1 g = t Xo R x {k)g 
for all k with Imk > and all g G L 2 (X ) which are in a neighborhood of dX . 
Since such g are dense in L 2 (Xo), this must in fact hold for all g G L 2 (Xq). 

Since (k 2 — Hcg)^ 1 = l Xo R x (k)i Xo for all k with Irak > and since both sides 
have meromorphic continuations to A -(A aXo ) (see [10] and Lemma [4. ip . they must in 
fact agree for all k G A (J (a 9Xo )- □ 

Lemma 3.4. Suppose (k 2 — iJcfj)^ 1 exists. For f G H 1 (dXo), let 

u=(k 2 -H cS )- 1 W(k)f. 

Then 

(3.7) (k 2 + A Xo )u = on X° , (9 n - iP 2 )u\ dXo = -P k f. 
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Proof. We first claim that there exists P G H (X ) such that 
(3.8) (d n -iP*)F\ 9x = -P k f. 

In fact, for < /i< 1 define N(h) : H s (dX ) -> H s ^(dX ) as follows: 

iV(/i) = i<J - A aXo )^ , iV^)" 1 = 0{h) : H s -\dX ) — > P S (<9X ) . 

Let x G Q)°([U; e)) be equal to 1 in a small neighbourhood of 0, with e chosen so that 
X ~ (— e, 0} x x <9X , near the boundary. We define (note that x < 0) 

[T{h)g\{x, y) = x(-x) exp(x(I - A dXo ) */h)g(y) , 
T(h) : H s (dX )->H s+l i(X ), s>0, 

so that 

T(h)g\ dXo = g , d n T(h)g\ dXo = N(h)g . 

For a fixed fc, P| = 0(1) : H 3 ^ 2 (dX ) — > H 1 ^ 2 (dX ), and hence, if /i is small enough, 
we have the following inverse 

(iV(/i) -iP^ 2 )^ 1 = iV(/i)- 1 (/-?P fc 2 iV(/ i )- 1 )- 1 : H^dXo) — >P§(<9X ). 

Using this and the mapping properties of T(h) we construct 

P ^ -T(h)(N(h) - iPlr'Pkf G H 2 (X ) , 

which satisfies f !3.8j) . 
We now set 

u = jP _( A; 2_ i/cfr) -l (A; 2 + Axo)jP5 

and observe that u satisfies the equations (13.71) . It remains to show that v = u. 
To see that we let h G C°°(Xo), and apply Green's formula to compute 

{{k 2 - H cS y l (k 2 + A Xo )P, h) = ((k 2 + A Xo )P, ((/c 2 - Peff)" 1 )^) 

= / (a n F((fc2 - H cS )~i)*h - Fd n ((k 2 - H cS )-i)*h 

JdXo V 

+ (F,(k 2 + A Xo )((k 2 -H cS y 1 yh) 

= [ (d n F((k 2 - Peff)" 1 )^ - Fd n ((k 2 - H cS )-i)*h 

JdX V 

+ {F,h). 

Now we use that d n F\g Xo = iP 2 (F\Q Xo ) — P k f, and that 

w G P 2 (X ) n Domain((& 2 - P cff )*) =>■ d n w\ dXo +i(P k 2 )*Rv = . 
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Thus we have 

((k 2 - H cS )-\k 2 + A Xo )F,h) 

= [ ((iP 2 (F\ dXo ) - P k f)((k* - H eS )-i)*h - iF((P*)*(k* - HaWk) + (F, h) 

JdX v 7 

= -/ P k f((k* - H cS )-i)*h + (F,h) = [ (-(k 2 -H eS )- 1 S Xo P k f + F)h, 

JdX Jx 

where the last expression follows from the definition of Sg Xo . Since this holds for all 

h e C°°(X ), 

v = F - (k 2 - H eS )-\k 2 + A Xo )F = (k 2 - H eS )~ l 6 Xo P k f = u , 

proving the lemma. □ 

We can now state and prove the main result of this section. It provides a justifica- 
tion of (JT5} and (TO) . 

Proposition 3.5. Let W be given by (13. ip . Then the XX' entry of the scattering 
matrix defined in £}H ^ s given by 

(3.9) S x ,y(k) = (Sf(k)ip x , <^a')l 2 (9x ) , 

where 

s f (k) = - (i - 2tw(kY(k 2 - H cE y l w{k)) , 

and H e s is defined in Lemma \3.B, 

Proof. We use Lemma I3~4l to express the action of {k 2 — H c s)~ 1 W(k). Suppose v\ = 
{k 2 — H c fi)~ l W{k)Lp\. Let U C X be a neighborhood of dX . On U we may use 
coordinates (x,y), with y EY. Since v\ lies in the null space of —A Xo — k 2 , we have 
that 

X' 

The boundary conditions (13. 7h applied to v x at dX mean that 

^ikx'iax' - bx>)ipx< - i ^2 kx '( ax ' + h'fyx 1 = -Pk^Px- 

X' X' 

Then bx = l/(2iy/k~\) and by = if A' ^ A. Thus v\ is the restriction to Xq of 
-i$ A /2, where $ A is determined by (12. ip and (I2.2p : 

(3.1UJ 3 > Al(o,oo)xy= e -= + 7 ^A'Ae 



. ^a ^ \kx' 
Therefore 

2 2 



= y/kxa\«px> - r<£>A = -77 ^ Syx'Px' + V?A 
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which proves the proposition. □ 

The equation (13.91) is valid for all real values of k (that is, k on the boundary of the 
physical space) with k 2 > a 2 , erf,, since the matrix coming from the right hand side 
is unitary and hence the singularities of (Sf(k)<p\>, (p\)L 2 (dx ) resulting from poles of 
(if c ff — k 2 )~ l are removable. 

4. Accuracy of Approximations 

Here we investigate the accuracy of the approximations made to use (I1.5P in nu- 
merical computations. Set 

nf °/ = l [OA] (-A dXo )f for / G L 2 (dX ), 

11^= l [0 , N] (H in )g iorgeL 2 (X ). 

In parallel with this, we introduce 

Woo,oo(^) *= W(k) , W^k) = WUf , 

W NjA (k) = I§W(k)n° x = I^WocaW , 

and 

= H cS , H N , A = H m - iW^W^, N G R U {oo} . 

Although Wn,a, H^oo.a depend on k, for simplicity we generally omit this in our 
notation. Note that H e Q, and -ffjv.A also depend on k. A quadratic form 

argument (see Lemmas 13.21 and 14.41) . using the form domain H 1 (Xq), shows that if 
u is in the domain of A , then d n u — iP 2 Tl a ^ x ° Ru = 0. However, for N < oo the 
domain of i?jv,A is the set of elements of H 2 (Xq) which satisfy the Neumann boundary 
condition, d n u = 0. 

Likewise, we define the approximations of the (full) scattering matrix obtained by 
using the approximation H^,a of H e s by S/,iv,A : 

(4.1) S f , N , A (k) = ~(l- 2zW N , A (kY(k 2 - H^y'WNAk)) • 

In order to bound the error in these approximations, we shall first see how close 
S/,oo,a is to n^ S , / iOOjOO , and then study the difference 

K (SW - s f , N , A ) n*f ° . 

4.1. Projection on 8Xq. We first analyze the approximation with a finite A and 
N = oo. The spectral cutoff for the boundary Laplacian, A has to be taken large 
enough to guarantee that Im k\ > for a 2 > A. The errors then come from evanescent 
modes and can be estimated using exponential decay. We present the results in two 
lemmas. 

Recall that H e s = H e s(k) is defined for k G A ct (a sy ). 
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Lemma 4.1. Let A CT (A aXo ) be the Riemann surface, given in (12.31) . to which the re- 
solvent of —Ax, (—A x — A; 2 ) -1 , has a meromorphic continuation (see [TU| Sect.6.7]j. 
Then the operators {k 2 — H^y 1 and (k 2 — -ffoo,A) _1 are meromorphic on A a ^ gXQ y If 
k 2 — H cf[ = k 2 — H eS (k) is invertible, so is k 2 — H^ ^k) for A > A (A;o) sufficiently 
large, and 

\\(k 2 - H^y 1 - [k 2 - < Ch." 1 ' 2 , A > Ao(fc) . 

Moreover, for k restricted to a compact set K C A CT (A aA - ) on on which k 2 — H c q is 
invertible, A and C can be chosen independently of k. 

Proof. Recall that U is a neighborhood of dX which we may identify with (— e, 0] x x Y 
with g\ U = (dx) 2 + gy. Choose Xi £ C°°(X), i — 1,2, so that each x% nas support 
in U, Xi = 1 in a smaller neighbourhood of the boundary, and 

X1X2 = Xi , supp x' 2 n supp xi = ■ 
Set R^ e (k) to be the operator on L 2 ((— 00, 0] x Y) defined by the Schwartz kernel 

R A , e (k)(x,y;x',y') d ^ £ -L( e **l»-*'l + (1 - nf °)e ik ^)cp x (y)cp x (y'). 

x x 

Note that R Afi (k) is a meromorphic function of k G A ct (a 9Xo ) since k\ is holomorphic 
on A CT ( AaA . o ). Let Roo, e (k) be the operator with Schwartz kernel given by 

and set 



2ik x 



E A (k) = (1 - xi)(A; 2 - H^y 1 + X2i?A, e (fc)xi , AgIU {00} . 

Then, for the same values of A, E A v satisfies the boundary conditions of if^A, that 
is 

(d n -iP*n 9 A Xo )E A v\ 9Xo =o, 

and is meromorphic on ). Moreover, 

{k 2 + A Xo )E A (k) = 1- [A Xo ,xi](k 2 - H^Y 1 + [A Xo , X 2]RA,e(k)xi 
= I + K A (k) 

where K A (k) is a compact operator. Moreover, K A (k) is a meromorphic function of k 
in A ct (a SYo ) with finite-rank poles. When k G iM + , ||i^A(^)|| — > as fc 2 — > —00. Thus 
/ + K A (k) is invertible for k G —k 2 ^> 0, and by analytic Fredholm theory (see 
for instance [T3l §2.4]) we have that 

{k 2 - H^k))- 1 = E A (k)(I + K A {k))~ l 

for k in the physical space, and it has a meromorphic continuation to A ct (a SXo ). 
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Now 



\E A (k) - E^k)^^ < CW^I" 1 . 



al>A 



For k in a compact set of A^a^ ) and o~a > A > A (fc), sufficiently large, we have 
Im k\ > 0, and since 

x £ supp x'2 ) x ' e supp Xi =>• I a; — x'\, \x + x'\ > e > , 

we have 

\\K A (k) -KniftW^v < Cmail^e- £ » Im ^ 2 . 



This constant is independent of k. Thus, if A is big enough, / + K A (k) — K 00 (k) is 
invertible with small norm, and 

|| (k 2 - H^y 1 - (k 2 - iWr 1 !! < CAT 1 ' 2 

for A sufficiently large (depending on k or K, eo and || (k 2 —H c s) The constant can 
be chosen independently of k on a fixed compact set K where k 2 —H C R is invertible. □ 

Remark. Using this Lemma and the definition ( 14. ip of 5/ )CO) a, we can see that for 
A e M+ U {00}, 

(4.2) p- l S fj00jA (k)P k = -(I - 2tP- 1 W O0 , A (k) t (k 2 - /Wr^A^) 

has a meromorphic continuation to A cr (A eXo ). The conjugation by Pt is necessary 
because while Pj 1 is a well-defined operator for k £ A ct (a SYo ), Pfc is not. Thus the 
operators W / o,a(^)-^' an( ^ PjT ^q^a^)* ar e well-defined on A CT (A aXo ), while in general 
W / oo,a(^) and W / ^ )A (/c) are not. The existence of the meromorphic continuation of 
(14. 2 p means that (v / ^v/v / ^a)(>S'/,oo,a(^) ( ^a'; fx) has a meromorphic continuation to 

A <x(A eXo )- 

Lemma 4.2. Fkt A < 00 and & so that k 2 — H c g is invertible and lmk\ > if 
a\ > A . Suppose / £ L 2 (<9X ) safe/^es n^ / = f for A > A . TTien, /or A > A 
swc/i t/iat fc 2 — -£?oo,a invertible, we have for some e' > 

||P fc - 1 ng ro (5 / (*) - S / , 00 , A (fc))P fe /|| L2(a xo) < Cmax(|A; A |exp(-e'Imfc A )) ||/|| l2{9 a„). 



1>A 



In particular, by Lemma\JJj this holds for all A sufficiently large depending on k. We 
note that the constants C and e' can be chosen independently of k if k is restricted to 
a fixed compact set K on which both of k 2 — H e g and k 2 — A are invertible and for 
which Imk\ > when a\ > A . 

We note that since n^ x ° / = /, the H 3 ^ 2 norm of / is bounded by a A-dependent 
multiple of the L 2 norm of /. 
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Proof. For A G [A , oo) U {oo}, set u\ = ik 2 — H^ ^^W^^Pkf ■ That means that 
ma satisfies 

(k 2 + A Xo )u A = on X° , (d n u A - iP 2 U a A x °Ru A ) = -P 2 f. 

Choose x £ e/2, 0]) to be one in a neighborhood of 0. Since U C X is 

a neighborhood of dX which can be identified with (— e, 0]^ x Y y , we can consider 
X — x( x ) to be defined on X by extending it to be outside of U. For g G L 2 (X), 
define U A x ° X g G L 2 (£7) C L 2 (X ) via 



Then 

(4.3) u A = (1 - x)«oc + nf O x«oo + (A: 2 - H^Y^k 2 + A Xo )(l - lif °) XUao 

since the function on the right satisfies the same boundary conditions as u A and is in 
the null space of k 2 + A Xo ■ 
Note that by using U 9 A X ° f = f 

(4.4) Uoo \u= (W*** + b x e ik ^)cp x + ]T b x e**<p x 

for some constants a\, b\, so that, using orthonormality of <f X s, 
hocWlz > Wu^u \\ 2 L2 > f \bxe ik " x \ 2 dx 

(4.5) = / \he^ x \ 2 dx= 



2Imfc> 



Also, 

(k 2 + A Xo )(i - nf «)x«oo = [dl - nf 



OO ) 



where x has the same properties as x an d XX = X- Our argument below takes 
advantage of the fact that the support of [<9 2 ,x] is contained in [e/2,0], while the 
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expansion (14.41) is valid for x in (— e, 0]. Hence, 

\\{k 2 + A Xo )(i - nf °) XWoo || 2 = ll^xK 1 - nf °)xuoo\\ 2 



< C( e " 4 ) / E (kx) 2 \bxe ik " x \ 2 dx 



; Im fc^ 



2Imfc 



e Im fc^ 



af>A 

Thus ffl~5l) gives 

||(fc 2 + A Xo )(l - ^xuooll < C(e- 2 )|| Woo |U 2 max(\k x \e-^/ 2 ) 



Using (14. 3p . the estimate 

\\(k 2 - H^^gWm < (1 + \k\)\\(k 2 - H^^gWv 
and the previous lemma, we obtain 



n d A fR( Uoo -u A )\\ LHdx) 

C(e~ 2 + |A;|)||(A; 2 - H^)-^ IKIb max (|fc A |e- eImfc ^ 2 ) 



< 



*1>A 



Thus far each constant C can be chosen independent of k, though of course ||woo|| 
depends on k in a continuous fashion on compact sets on which k 2 — H e R is invertible. 
Note that P^PkU 9 / = U 9 A X ° is a bounded operator. Thus using the expression for 
Sf, S/,oo,A and the previous lemma finishes the proof. □ 

4.2. The cut-off in the interior. We now turn our attention to the error introduced 
by using LT^. Throughout this section we assume that A < oo. 
Our results will use the following standard 

Lemma 4.3. Suppose X is a compact Riemannian manifold without boundary and 
X £ C£°(M) is equal to 1 in a neighbourhood ofO. Suppose that Y C X is a smooth 
embedded submanifold of codimension one. Then 

11(1 - x(-h 2 Ax))u\ Y \\lhy) < CVh\\u\\ Hl ffi . 

IfveH 2 (X\Y)D H\X) then 

||(1 - x(-h 2 Ax))v\ Y \\l*(y) < Ch\\v\\ H2{ x\Y) ■ 
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Proof. Both statements in the lemma are local. In fact, if P is another elliptic sec- 
ond order operator on X then for some constant Cp the calculus of semiclassical 
pseudodifferential operators (see for instance [U Appendix E]) shows that 

(1 - xhh 2 A))(l - X (-h 2 C P (P + C P ))) = (1 - x(~h 2 A)) + H - k ^ Hk (h N ) , 

for all iV and k. Hence we can use any other second order elliptic operator and that 
property is invariant under changes of coordinates. 

It follows that we can assume that X = M. n and Y = {x± = 0}, R™ 3 x = (xi,x') 
(the compactness is irrelevant for the local statement). 

Denoting the Fourier transform by T we write 

(4.6) ((i - xC-^A^Ky) (O = / (i-x(^ 2 |£| 2 )H£i,£'Ki- 

Jr 

Hence, by the Cauchy-Schwartz inequality, 

ii(i - x(-h 2 Ax))u\ Y ni 2(y) < c [ h)(i - x(/^i 2 ))koi 2 (i + \en , 



where 



F(e',/i) d = / (l-x(^ 2 |e| 2 ))(l + le| 2 )- 1 ^i 
Jr 

< I (i + lai 2 )- 1 ^ + %i>ca(o / (\a 2 + i^rr 1 ^ 

J\ti\>c/h Jr 

< Ch. 

This proves the first part of the lemma. 

For the second part, we can assume that suppt> C {x G R n : |x| < R} as we can 
localize to a compact set. We then write 

r R 

(4.7) v(0 = / (e- ix ^^ x ^v( Xl , O + e^F^vi-xu £')) dx x . 
Jo 

Since v G H 1 (W 1 ), J r x ' h ^^v(0,C,') G L 2 (M n_1 ) is well defined and an hence we can 
integrate by parts to obtain 

Since v G i/ 2 (lR^), c^i^Oi, £') is well defined in L 2 (M n_1 ). We now use the following 
decomposition: 

(i - X (h 2 e))v = vi+h, mo = ht>c/ h (m - x(h 2 e))m , 
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noting that > c/h on the support of w 2 (£). We first estimate the contribution of 
t>2 as in the proof of the first part of the lemma: 

IMy \\h<y)<C f G(g,h)\m\ 2 0- + \e\ 2 ) 2 dt 



where 



/i) = f / (i-x(h 2 \em+w\ 2 r 2 dti 

J\b\<c/h 

p2c/h 

< %i> cA (0 / (i + ief)- 2 d6 

JO 

< Ch\ 

which is a better estimate than needed. 

To estimate the contribution of Vi we use ( 14.71) : 



\vi\ Y \\h( Y ) < C r\\ v W 2 H 2 ( x\y) [ / TI^i ) < C R h 2 \\v 

\J\£i\>l/h ?1 / 



2 

2||„.l|2 

H 2 (x\y) ' 



'|ft|>i/fc £i 

which completes the proof. □ 

Like H c s, Hn,a = -£^v,a(&) is a well-defined operator for k £ A -(a SXo ). 

Lemma 4.4. Fzx A < oo, and suppose that k 2 — ^ is invertible, k £ A tT (A aXo )- 
Then, for N sufficiently large, k 2 — H^,a is invertible, and 

\\(k 2 - H^r 1 - (k 2 - ^,A)^ 1 |k-i(x„)- > ^(x ) < CN- X '\ 

The constant C can be chosen uniformly for k in a compact set K C ^a{A aXo ) on 
which k 2 — i?oo,A is invertible. 

Proof. As in £J3] we will use quadratic forms to reinterpret our operators. Thus, for 
N eR + , E eC, set 



?OO)A 0, £)(«, u) = / VuWv - i \ W^ A uW^ A v — E uv 

VuW-i / P k U d A Xo RuP*U d A Xo Rv - E 

o JdX 



and 



/ uv 

lx 



q N>A (k, E)(u, v) = VuVv - i / W l NA uW* NA v - E 

JX JdXo 

= [ / P k U d A Xo RU%uP^U d A Xo Hn.^v - E [ 

JXn JdXn 
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Here we take both form domains to be H 1 (X ). The quadratic forms q^^kjE) 
and qN,\{k, E) are associated to operators -ffoo,A — E and H x ,a — E respectively. We 
expand the difference of the quadratic forms as follows 

g^A (&,£)(«,«) - q N) A(k,E)(u,u) 

= [ {\p k u 9X °R U \ 2 - \p k u 9X °iffl>| 2 ) 

JdX 

= [ ((P k U dXo R U )P k U d A Xo R(I - W*)u + P k U 9X °R(I - li^ U P k U dXo RU%u) . 
Jax Q v ' 

We have the following estimates: 

||nf °r(i - n™)u\\ m{dXo) < c £jA \\r(i - n%)u\\ LH9Xo) < c^n-^Mh^ 

(4-8) Qx 

||I1 A °Ru\\ H e( dXo ) < C^A||_Ru|| i 2( eX() ) < C^AlMltfipfo) ■ 

To obtain the first, we apply Lemma 14.31 to X = f X U Xq where the metric on X 
is obtained by reflecting the metric on Xq through Y = 8Xq. Since the metric has 
product structure near Y this means that 

(here ev refers to even functions) and the action of the Neumann Laplacian on X is 
the same as the action of on even functions. Applying Lemma l4~3l with h = l/\/~N 
gives (HSJ). 

Applying ( 14. 8ft to estimate the difference of the quadratic forms we obtain, for 
E < 0, 

\ qoo>A (k, E)(u, u) - q N>A (k, E)(u, u)\ < C K (k)N- l ' 4 \\u\\ 2 Hi 

< C K {k)N- l ' i Re{q ooA {k,E){u,u)) . 

The constant depends continuously on k. Here we use the fact that Im k\ > for all 
but finitely many A, ensuring that Req^^k, E)(u, u) bounds IMI#i(x ) f rom above 
for £<0. Thus, by [BJ Theorem 3.4], 

II (-#«>, a - EY 1 - (H NA - ^"l^CXoHL^Xo) < CN" 1 ^ 
for N sufficiently large depending on E, k, and A. This dependence on k is continuous 
on regions where k 2 — f/oo.A is invertible. To extend this to other values of E (in 
particular, E = k 2 ), we use 

(4.9) 

(A - z)- 1 = {/-(/+ ( z - Zo )(B - z)- 1 ) (z - zo) ((A - zo)- 1 - (B - zo)" 1 ) J" 1 

x (/ + (z - z )(B - z)~ l ) (A-Zo)' 1 . 
Consequently, if k 2 — i^oo,A is invertible, so is k 2 — Hn,a f° r sufficiently large N, with 
II (k 2 - tfocA)" 1 - (k 2 - H n> a)- 1 \\l*{x )^{x ) < CN- 1 / 4 
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Here the constant will depend on k and A, as will the lower bound on the N for which 
this holds. These can be chosen uniformly if k is restricted to lie in K. 

Now we show that there is a similar bound from H~ 1 (Xq) to i7 1 (Xo). We choose 
E so that Re(q OCJ: \(k, E)(u,u)) > Co\\u\\jji(x ) f° r some c > and all u G H l (X ). 
Suppose w G L 2 (X ) and set it = (i^oo.A — E)~ 1 w, un = {Hn,k — E)~ lr w. Then 

c 1| < Re (qoo,k{k, E){u, u)) = Re(w, u) 

so that 

colMl^i < IMI-ff- 1 - 

This shows that we can (uniquely) continuously extend (Hoo^ — E)' 1 to be a bounded 
operator from H~ 1 (X ) to H 1 (X ) when £ < (the duality argument shows that 
we can extend the operator to the dual of H 1 (Xq) and H~ 1 (Xq) is contained in that 
dual as the space of elements of H~ l (X) supported in X ). The resolvent equation 
extends this to other values of E. Likewise, 

c \\u - u N \\ 2 H i < Re (qoo tA (k, E)(u - u N ,u - u N )) 

= Re(q 00j A(k, E)(u,u - u N ) - q^^k, E)(u N ,u - u N )) 
= Re ((w, u — un) — (w, u — un) 

+ q N ,A(k, E)(u N ,u - u N ) - qoo s A(k,E)(u N ,u- u N )) 

< CA^~ 1//4 ||MAr||^i \\u — Ujy\\ H i 

< CiV" 1/4 (||n|| // i + ||u - u N \\ H i)\\u - u N \\ H i. 
We allow the constant C to change from line to line. This implies that 

\\u ~ u N \\ H i < CiV _1//4 (||M|| f/ i + ||m - u N \\ H i), 

which then means that for sufficiently large TV 

\\u-u N \\ H i < l7A^- 1/4 ||m|| h i. 

Using (14. 9 p this can be extended to other values of E. Again, these constants can be 
chosen uniformly for k G K. □ 

Lemma 4.5. Fix A < oo and k so that k 2 — a is invertible and lmk\ > if 
a\ > A. Suppose f G L 2 (dX ) satisfies U 9 / f = f. Then, for N so that k 2 - H N>A 
is invertible, there is a constant C depending on A and k so that 

\K x "P^\s Loo Ak) - s LN Ak))P k f\\ < CN-^WfW^y 

The constant C can be chosen independently of k, if k is restricted to a compact set 
K C A (7 (a (9x ) on which k 2 — ^ and k 2 — H N A are invertible. 
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Proof. Choosing N so that k 2 — iijv,A is invertible, set 

Uoo = {k 2 - Hoo^^Woo^Pkf and 
u N = (k 2 - H N ^y l W NA P k f. 

That is, satisfies 

(k 2 + A Xo ) Uoo = on X° 
dnUoo - iP^U dXo R Uoo = -P 2 f 
and un satisfies, for N < oo, 

(k 2 + A Xo + iW^aW^uj, = W N , A P k f on X° 
d n UN\ax = 0. 

Note our assumptions on / mean that the if 3//2 norm of / is bounded by a A dependent 
constant times the L 2 norm of /. 

We wish to understand IT]y'u 0O . Let ^f n be a real eigenfunction of the Neumann 
Laplacian on X , with —A Xo ^J n = T 2 ^ n . Suppose in addition that H^nlU 2 ^) = 1- 
Then 

(t 2 - /c 2 )(Moo, ^n)L^(X ) = ~ (A Xo ^ n Uoo ~ ^ n A Xo Uoo) dvo\ Xo 

^ndnUoo dvoly 

dX 

= [ y n (lP 2 n 9 A X «R Uoo - P 2 f) dvoly. 
JdXo 

-A Xo - k 2 )U^ Uoo = iW N , A W^ AUoo - W NA P k f. 



That is, 
Thus 



u N = n^noo - (k 2 - H^y 1 ((k 2 - H NA )U^ Uoo - W NA P k f) 
= U^ Uoo + iik 2 - H NA y 1 W NA {W t ooA - M^ iA ) Ubo 
(4.10) = IT^oo + i(k 2 - H NA )~ 1 W N , A P k n 9 A Xo R(I - n%) Uoo . 

The second part of Lemma 14.31 gives the following estimate: 

\\R(1 - n£) Uoo || L 2 (aX) < CN" 1 ' 2 llMooUfP^o) , 

and consequently, 

\\p k u 9 A Xo R(i - n^uJl < cn-^Wu^Wh^ 

with constant C depending continuously on k. 

Let g e L 2 (dX ), h e HV 2+ (X ). Then Rh e L 2 (dX ), and 

\(W NA g,h) Xo \ = \(g,P* k U d A Xo RU^h) dXo \ < \\g\\m d x Q )\\h\\ H1/2+ {Xo) . 
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That is, 

||W / at,a5 , IIh- 1 (^o) < I|Wjv,a<7||#-i/2-(x ) ^ C|blU 2 (x ) 

and the constant is independent of N and depends continuously on k. On the other 
hand, Lemma 14.41 shows that for TV sufficiently large 

|| (k 2 - H NA )~ ||if-i(x )-»Hi(x ) < C. 

Using these estimates in (14.101) . we find that 

Ikoo - u N \\ H i < CN~ 1/2 , 

implying the desired bound by restricting to 8Xq and using again the fact that 
P^PkU 9 / = U 9 A X ° is a bounded operator. □ 

5. Proofs of Theorems 

Our proof of the Theorem in £Q] wn l use the unitarity for k real not only of the 
finite-dimensional scattering matrix defined by ( 11. 71) . but also of the approximations 
of the scattering matrix obtained by introducing the projections IT^ and II 9 ^ . 

Lemma 5.1. Let k e R. Then S(k) defined by (T7p and U 9 k x ° S f>N)A {k)U d k x ° , for 
A, N G M+ U {oo} ; are unitary. 

Proof. That S(k) is unitary for k real is well known. It can be seen as follows. Recall 
that S(k) = U 9X ° Sf j00)00 (k) n|f o . 

We note that (P k )*<px = k\(p\ and k\ is real for o\ < k 2 and pure imaginary for 
a 2 > k 2 . Therefore 

(5.1) U 9x »W\k) = U 9X °W*(k) and (/ - 11**°)^) = -(I - n d *°)W*(k). 
Thus, we have 

(5.2) W{kY*U 9 A Xo W*{k) + W{k)Ii 9Xo W\k) = 2W{k)U 9 A Xo U 9 k Xo W t {k). 

Using this and the resolvent identity gives 

(k 2 - H^)- 1 - (k 2 - H*^)- 1 

(5.3) = %{k 2 - H^r^-U^W^U^W*^ - U^WU 9X <WlI£)(fc 2 - tf^)" 1 

= _2i(k 2 - H^-^wul^u^w'u^k 2 - H* NA y\ 
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Therefore, 
(5.4) 

= nf° (7 - 2iw t NA (k)(k 2 - H^y'w^Aik) 

x Uf° (/ + 2zW* NA (k)((k 2 - H^r'YiWUk))*)^ 
= IT^ (/ - 2iW t NA (k) [(k 2 - H^a)- 1 - {k 2 - H^y 1 } W NA (k) 

+ AW^ik 2 - HM^WN^l^W^^k)^ 2 - H^-'W^ik^U 9 ^ 
where we have used (15.11) . Applying the identity (15.31) we find that 

(5.5) u d J o s LN Ak)nf%u d J o s f , N Ak)u d J°r = I 

as desired. □ 

Theorem 2. Let X be a manifold with infinite cylindrical ends, and Sxx>{k), Sf t N^(k) 
be as defined via 112. h2.2\) and fl^.ip . Suppose k G Iy a ^ gx ) and A G R are such 
that k 2 — H c ff = k 2 — H c g(k) is invertible, and Imk\ > if a 2 > A . Then, for 
o~l, o~\i < A and A > A , 



\fk\' 



S xy (k) = (P^S f , NA (k)P k ip x ,, (fx) + 0{N~* + e~ A / c ) . 



We recall that k 2 — H e s is invertible if k is in the physical space with Imk > 
0, ImA^ > for all A, and that (k 2 — H c fi)~ l is meromorphic on A tT (A l9Xo ). 

Proof. The proof follows from writing 



S X y(k) = (p- 1 S f (k)P kV >y,<p X ) 



(5 ' 6) = (P^3 f ,NA^)Pk<PX', <Px) + (P^iSfik) - S ftO0A (h))P k <p X r, if X ) 

+ (P^ 1 (5'/ )0Oi a(A;) - S fyN) K(k))P k (p X ',(f\) , 

where we note the first equality follows from Proposition 13.51 Applying Lemma 14.21 
and Lemma 14.51 we obtain the theorem. □ 

We now prove Theorem [IJ 

Proof. If k 2 — H D g is invertible this is just Theorem [21 and hence it remains to prove 
that the estimate is valid for all fcfl, even if k 2 — H e g is not invertible. 

Using the unitarity proved in Lemma loTTj along with the fact that a 2 , < k 2 , o\ < k 2 , 
we see that each of the terms on the right hand side of (15. 6p is bounded for all 
k G M. Also, for N, A G M + U {oo} (\/ky/ \/k~x){Sf t N,A(k)ipx, ipy) has a meromorphic 
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extension to A cr (A ex ), as can be seen from the formula (14. ip and the fact that (k 2 — 
-£^v,a) _1 continues meromorphically to A a ^ gx y Hence 

k G (-ax, a x ) H (-ay, cry) , 

has a neighborhood in A a ^ gx ) on which (Sf n Afx', fx) is holomorphic, N,Ae M + U 
{oo}. 

We will now apply the maximum principle: (15. 6p and Lemmas 14.21 and 14.51 show 
that 

S\\> — (Sf j N,A<PX',<P\) , 

is bounded by C(N~^ +e~ A /°) on the boundary of the neighbourhood chosen above, 
since (k 2 — H^)^ 1 is bounded there. The theorem follows as the difference is holo- 
morphic. 

In other words, we have shown the theorem holds when k G K. is on the boundary 
of the physical space even if k is a pole of (k 2 — if e fj) _1 , as long as k 2 ^ af, a 2 ,. 

To finish the proof, consider what happens at a point k^ G M, k^ = a 2 > a 2 , . (The 
case a\ < a\, follows by symmetry.) If a\ = af,, then since \fk~x' /\/k~x — 1 (except for 
the removeable singularity at kx — 0), (Sf^,A(k)<fx', fx) has a meromorphic extension 
to a neighborhood of k in A CT (A aXo ). The boundedness at k , again obtained from 
unitarity, ensures that there exists a neighborhood of k on which (Sf ) N,A<Px , ,<P\) is 
holomorphic. Thus the previous argument using the maximum principle holds here 
as well. 

Now suppose a\, < a\ = k\, and set T\y(k) = (y/ky /y/k\)S\y(k). Then T\y is 
meromorphic in a neighborhood of k$. Using the unitarity of S(k) for k real, 

\Tx\>(k)\ 2 < l^vll^r 1 for k2 > ^A2, keR. 

Thus y/k\T\y(k) is bounded at k , and T X y must then also be bounded at k , since 
near ko it is a meromorphic function of k\. Therefore S\y(ko) = 0. Since we have 
in fact only used the unitarity of S(k) for k G 1R and the existence of a meromorphic 
extension, the same argument gives 

(Sf, N ,A( k ) ( Px'^x)\k=k a = 0; for a\, <a\ = k 2 ; iV,AeM + U {oo} . 
Thus the approximation is exact in this special case. □ 

6. An example 

In this section we consider the simplest one-dimensional example where things are 
explicitly computable and we are able to see the effects of the approximation 11^ 
explicitly. 

Let X = (— 7r, oo), with X = (— vr, 0] and X\ = [0, oo). We consider the operator 
—d 2 on X, with Neumann boundary conditions. Although strictly speaking this 
example does not fall in the class considered in the first part of the note (X has 
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Approximation with M-10 




Approximation for different values of M 




Figure 2. An illustration of the example in £j6j the top figure shows 
the real parts of the approximation and of the scattering matrix, and 
the lower one, the graphs of \SM 2 {k) — S(k)\ for different values of M. 

a boundary, {—it}), it is easy to see the arguments of the previous sections follow 
through, with dX replaced by Y = {0}. Because Y is a point, the full scattering 
matrix is a scalar, and is easily computed to be S(k) = e 2mk . 



For this example, 



71 



-1/2 



if n = 



(2/7r) 1/2 cos(nx) if n > 0. 

Since there is no sense in the cutoff n^ x ° for this problem, we use only one subscript 
on our approximations of W: 

M 

w M 2(k) d = VkJ2^(o)y n (x). 



n=0 



Similarly, we denote the approximation of S(k) thus obtained by SM 2 {k)- In the 
notation of the paper M = s/N ' . We denote by Wm 2 = W^pik) the M + 1 vector 
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Note that if ia l a 7^ — 1, 

(I + iady 1 — I — + ia t aY 1 aa t . 

Set D M 2 = D M 2(k) to be the M + lxM+1 matrix given by ((A; 2 — n 2 )5 nm ). We see 
that when k ^ Z so that DM^(k) is invertible, the approximation SM 2 {k) is given by 

S M2 (k) = -1 + 2zW^ 2 (£> M 2 + i^^)" 1 ^ 

= -1 + 2i(D^ /2 W^)*(/ + iD'^W^D-^WMtfT^^WM*. 

Set 5 M 2 = D~^ 2 /2 W M 2 and /3 M 2 = B t M2 B M 2. Then, for fc ^ Z, 

5 M 2(/c) = -1 + 2z^ /2 (J - z(l + ip M 2)- 1 B M 2B t M2 )B M 2 



(6.1) =-l + 2z(/3 M 2 
Now 

(6.2) /3 M2 
We note that 







M 2 




lim /5 M 2 = cot 7r/c; 

M^oo 



one can use this and (16. ip to see that 



lim S M 2(k) = e 2mk = S(k) 

when k 2 £ N . Using (jSTTJ) and (1Q]) . we see that for jfe e R \ Z and M > 

Ci/M < |S M2 (fc) - < C 2 /M 



for some positive constants C\, C2 depending on k. Since M = y/N, this shows that 
the estimates obtained in Lemma 14.51 and in the main theorem are optimal. 
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